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Abstract 

We examine the lowest order quantum corrections to the effective 
action arising from a quantized real scalar field in the Randall-Sundrum 
background spacetime. The leading term is the familiar vacuum, or 
Casimir, energy density. The next term represents an induced gravity 
term that can renormalize the 4-dimensional Newtonian gravitational 
constant. The calculations are performed for an arbitrary spacetime 
dimension. Two inequivalent boundary conditions, corresponding to 
twisted and untwisted field configurations, are considered. A careful 
discussion of the regularization and renormalization of the effective 
action is given, with the relevant counterterms found. It is shown that 
the requirement of self-consistency of the Randall-Sundrum solution is 
not simply a matter of minimizing the Casimir energy density. The 
massless, conformally coupled scalar field results are obtained as a 
special limiting case of our results. We clarify a number of differences 
with previous work. 
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1 Introduction 



The search for a unified field theory led from the very early days after 
Einstein's theory of general relativity to the suggestion that there could be 
extra spatial dimensions beyond the familiar three. The original idea of 
Kaluza Q was to incorporate electromagnetism into a 5-dimensional theory 
of gravity. In order to explain why the extra spatial dimension has not 
been observed, Klein |Q] suggested that it was compactified to a circle of 
extremely small radius, typically taken to be of the order of the Planck 
length. The generalization to more than one extra spatial dimension [^] 
allowed the incorporation of non-abelian gauge fields as part of the higher 
dimensional metric. These higher dimensional theories are usually referred 
to by the generic term Kaluza-Klein theories. 

In the early 1980's an immense amount of interest was generated in 
Kaluza-Klein theories by Witten's Q observation that seven extra spatial di- 
mensions was the minimum required for incorporation of an SU{3) x SU (2) x 
U{1) gauge symmetry, and Nahm's proof ^ that this was the maximum al- 
lowed by supergravity. This focused attention on 11-dimensional theories 
with various compactifications. (See Ref. Q for a review.) Further inter- 
est in higher dimensional theories was, and continues to be, generated by 
superstring theory. 

A new angle on Kaluza-Klein theories was pointed out recently In 
contrast to the standard belief that extra dimensions must be associated 
with extremely small length scales, it was noted that the extra dimensions 
could be much bigger than the Planck scale and still remain as yet unde- 
tected. This has the effect of bringing the fundamental Planck scale much 
closer to the electroweak scale, and thus offers a possible solution to the 
gauge hierarchy problem. The basic idea behind this is to have a (4 + n)- 
dimensional spacetime, with the n-dimensional space having a volume Vn- 
The relationship between the 4-dimensional Planck scale Mpi and the fun- 
damental mass scale M of the higher dimensional space is Mpi ~ V^M"+^. 
For large volumes Vn, we can have M of the order of the electroweak scale 
(~ 1 TeV) and still have Mpi ~ 10^^ GeV. Unfortunately this scenario 
does not provide a completely satisfactory solution to the gauge hierarchy 
problem since it trades a large ratio between the Planck scale and the elec- 
troweak scale for a large ratio between the compactification scale V~^^^ and 
the electroweak scale. 

Randall and Sundrum |^] suggested a model that overcomes the problem 
just described. This model is based on a 5-dimensional spacetime with the 
extra spatial dimension having an orbifold compactification. Two 3-branes 
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with opposite tensions sit at the orbifold fixed points. The Hne element is 

ds"^ = e-'^'^'^^Si^udx^'dx" - r'^d(t? (1.1) 

with x^^ the usual 4-dimensional coordinates, \(j)\ < vr with the points (x^, (p) 
and (x^, —(f)) identified. The factor of e~'^^^^'^^ in the line element is often 
referred to as the warp factor, since it results in the spacetime not being the 
direct product of 4-dimensional spacetime and the extra dimension. The 3- 
branes sit at <^ = and (p = ir. /c is a constant of the order of the Planck scale 
(the natural scale for the theory), and r is an arbitrary constant associated 
with the size of the extra dimension. 

The interesting feature of the Randall- Sundrum model is that it can 
generate a TeV mass scale from the Planck scale in the higher dimensional 
theory. A field with a mass tjiq on the = vr brane will have a physical mass 
of m ~ e~'^^^'mo. By taking kr ~ 12, and mo — 10^^ GeV, we end up with 
m ~ 1 TeV. In the original version of the Randall-Sundrum model all of 
the standard model particles are supposed to be confined on the brane, with 
only gravity in the bulk (5-dimensional) spacetime. Alternatives to confining 
particles to the brane, and allowing fields to live in the bulk spacetime have 



been investigated. See for example M, 10, 11, 12, 13, 14, 15, 16, 17, 18 



It is the possibility of bulk matter fields that we are interested in this paper. 

The radius r is assumed to be the vacuum expectation value of a scalar 
field, called the radion. In the original Randall-Sundrum model the value 
of r is undetermined. In order to make the theory physically acceptable a 
mechanism for determining r must be found. In Ref. a classical method 
for stabilizing r was suggested. 

The similar problem of fixing the size of the extra dimensions was found 
in the older Kaluza-Klein theories that were based on direct products of 
Minkowski spacetime and a homogeneous space, often chosen to be a sphere. 
It was necessary to have the extra dimensions of the order of the Planck scale, 
but there was nothing in the classical theory that determined this. It was 
recognized by Candelas and Weinberg [21| that quantum effects from matter 



fields or gravity could fix the size of the extra dimensions in a natural way. It 
therefore is sensible to ask if a similar procedure can fix the value of r in the 
Randall-Sundrum model. This problem has received some recent attention 



[p2| , p^ , 24] with some differing conclusions. The main purpose of the present 
paper is to provide details of what we believe to be the correct and complete 
calculation of the vacuum energy density for bulk scalar fields, and to calcu- 



late the induced gravity term in the effective action |25| since it is known 



that this can play an important role in self-consistent solutions and can give 
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different conclusions from merely considering the vacuum energy density. 
We also clarify the two inequivalent field configurations which can arise as a 
result of the boundary conditions that can be imposed. Finally, we analyze 
in detail the requirement that the Randall-Sundrum spacetime be a solution 
to the quantum corrected Einstein equations. It is shown that this is not 
just the condition that the Casimir energy density have a minimum. 



2 Quantum corrections to the classical theory 

In order to calculate the quantum corrections to the classical theory we will 
follow the procedure used in Ref. |2^. This consists of first expanding the 
higher dimensional field in terms of a complete set of modes appropriate to 
the extra dimensions. The dependence on the extra dimensions can then be 
integrated out of the action functional to leave an equivalent 4-dimensional 
theory, but one with an infinite number of fields. The masses of the fields 
are usually quantized in some way. This is the procedure used in our com- 
putation of the vacuum energy density [^3| for the minimally coupled scalar 
field. 

If we wish to calculate the induced gravity contribution to the effec- 
tive action we must generalize the metric (1.1) to allow curvature in the 



4-dimensional submanifold representing the branes. We will take 

ds^ = e-^''^y^g^^{x)dx^'dx'' - dy'^ (2.1) 
where y is related to the angle 4> in (|l.lD by 

y = r(t) (2.2) 

and the function cr(y) is given by 

a{y) = k\y\ . (2.3) 

We will allow the submanifold specified by y = to have an arbitrary 
dimension D. (The Randall-Sundrum model is then specified by D = 4 and 

The Einstein-Hilbert gravitational action is 

Sg = (167rG)-^ j d^xdy \g\^''^ {R - 2A) (2.4) 

where we use a caret to denote a (Z) + l)-dimensional object. (So gjxy is 
the [D + l)-dimensional metric found from ( |2.lD , g is its determinant, and 
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R is the scalar curvature constructed from it.) There are also terms in 
the total action functional coming from the two branes at y = 0, vrr which 
represent the hidden and visible branes respectively: 

Sbrane = - j d"" xdy \g\^/^ {VJ{y - Trr) + VhS{y)} . (2.5) 

We will return to this part of the action in Sec. ^. 

In order that the metric in (^]^) with gfj,,y = rj^^ satisfy the field equations 



coming from varying ( p.4[ ) and ( p.5| ) with (T{y) given in (2.3) along with the 
Z2 identification (x^,y) ~ {x^,—y), it follows that 

a" = 2k [6{y) - 5{y - vrr)] , (2.6) 

Vn = -V. = (2.7) 

A = -]^D{D - l)k'^ . (2.8) 

These reduce to the results of @ when D = 4. 

For later use we record the curvature scalar for ( p.ll) : 

R = e^''R + D[2a" -{D + l){a'f] (2.9) 
= e^''R + AkD[5{y)-5{y-TTr)]- D{D + l)k'^ . (2.10) 

We will consider a single real scalar field ^ of mass m in the spacetime 

S=\ j d^xdy |^A^5^$5.$ _ ^2^2 _ ^^^aj ^2.11) 

We allow a general non-minimal coupling to the curvature with dimension- 
less coupling constant ^. The minimally coupled scalar field has ^ = 0. 
The conformally coupled scalar field has ^ = (Z? — 1)/(4L'). We will keep ^ 
arbitrary here. 

By considering the {D + l)-dimensional field equation for following 
the analysis of Ref. 1^], leads us to define fn{y) as a solution to 

+ACDke-^''[6{y) - 5{y - 7rr)]/„ = mlfn , (2.12) 

normalized by 

+7rr 

j dye^^-^^''fn{y)fn'{y)=Snn' ■ (2-13) 
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We have defined 



m 



m 



D{D + l)k'^i. 



(2.14) 



We first obtain the boundary conditions satisfied by fn{y) by integrating 
(|2.12| ) first about y = and then about y = vrr. We assume that fn{y) is 
continuous at these two points. This gives the two resuhs 




dUy) 



dy 



dfniy) 



y=e 



dy 




dfn{y) 



dy 



dfniy) 



y=7rr—-e 



dy 



4^kDUiO) , (2.15) 
4^kDfn{7rr), (2.16) 



y=—7Tr—e 



(In the second result we have used the y ~ —y identification.) 

Next we solve ( |2.12| ) for y > {y ^ vrr) and y < {y ^ —vrr). Call the 
solutions in these two regions f^{y). It is easy to show that 



/n(2/) 



±ky 

k 



^ ±ky 
k 



)} 



(2.17) 



where and are constants of integration. The order of the Bessel 
functions v is given by 



L»2 



(2.18) 



At this stage we must choose how fniy) is related to fniy) by the Z2 
identification on 5^. The situation is like that used some time ago |26] in 
considering fields on non-simply connected spacetimes. For a real scalar field 
there are only two possibilities. We have either 



fni-y) = fniy) , 
corresponding to an untwisted field, or else 

fn{-y) = -fn{y) , 



(2.19) 



(2.20) 



corresponding to a twisted field. The first choice corresponds to the bound- 
ary conditions adopted in R ef. || ] and Ref. |22]. We examine each possibility. 
For the untwisted field ( ^.19 ) leads to 



A: 



in which case we can write 



fn{y) 



e 2 



k\y\ 



^An 



+ BnY,, 



(2.211 



(2.22) 
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for either positive or negative values of y (but y ^ 0, irr), as given in Ref. [^. 
Applying the boundary conditions ( p.l5| ) and (|2.16| ) results in the eigenvalues 
nin determined by solutions to the transcendental equation 

F„O=0, (2.23) 

with 

Fu{z) = - [y,{az)j,{z) - j,{az)y,{z)] , (2.24) 
vr 

where ju{z) and yu{z) are expressed in terms of Bessel functions by 

Uz) = ^Dil-AOJu{z)+zJUz) , (2.25) 

Vuiz) = ^Dil - AOY,{z) + zYliz) . (2.26) 

We have written 

a = e-^^' . (2.27) 

For = and D = A these expressions reduce to the results of Ref. |^] . 
For the twisted field (|2.2[l| ) leads to 

A- = -At , B- = -B+ (2.28) 



instead of ( 2.21 ). This time we have 

UV) = '^l^l [Ar^J. (^e^l) + B.y. (^e'^l^l) } (2.29) 

as the solution for ?/ > (y 7^ vrr) and y < {y ^ — vrr). Because fn{y) for 
the twisted field is an odd function of y, using the y ~ —y identification, we 
must have 

/n(0) = , fn{7rr) = . (2.30) 

These last two results are consistent with the general boundary conditions 
(|1|) and (|1|) since the derivative of an odd function is even, resulting in 
the left hand sides of the two equations vanishing. The boundary conditions 
( |2.3[1| ) lead to a much simpler analysis, with the eigenvalues m„ given as 
solutions to 

(^) = (2-31) 

with 

F^{z) = - [Mz)Y^{az) - Maz)Y^{z)] . (2.32) 

TT 
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Thus there can be no conflict with the results of Ref. |23] that used untwisted 
boundary conditions. As we discuss in Appendix B, there are some slight 
differences in the technicalities between the twisted and untwisted scalars. 

We now expand the (D + l)-dimensional scalar field ^{x,y) in terms of 
the fniy) as 

^x,y)=J2Mx)fn{y) ■ (2.33) 



This expansion can be substituted into ( p.ll| ) and the integration performed 
over y. The result is 

S = -\Y.I d''x\g\'/\n{^ +ml + (R)ipn . (2.34) 

This shows that we may view the theory as one in D dimensions, with an 
infinite number of fields whose masses are quantized according to ( p.23 ). 



At this stage we can work out the one-loop effective action in the usual 
way: simply compute the contribution from the v}^ mode, and then sum 
over all modes. Because we are after the vacuum energy, as well as the 
induced gravity term, it is advantageous to use a heat kernel method, and 
we follow Ref. |25] here. (See |2?] for a review.) 

The one-loop effective action for the n*'' mode is 

oo 

r« =-'-j d^:E|5|i/2 j ^ e-^^<K(s; x, x) (2.35) 



where 

d 

i—K(s;x,x') = (a+^R)K(s;x,x') , (2.36) 

OS 

K{0;x,x') = 6{x,x') . (2.37) 

The result for r^^'' is infinite and requires regularization. We choose dimen- 
sional regularization here with D viewed as D + e with e — > understood. 
We now make use of the asymptotic expansion (see Ref. Q for example) 

oo 

K{s; X, x) ~ it{4ms)-^/^ Yj^is)''ak , (2.38) 

fc=0 

where 

ao = 1, (2.39) 
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Higher order terms in ( p.38| ) involve higher powers of the curvature tensor; 
so the oo and ai terms contain the leading order terms. 
The two leading terms in F^^^ are found to be 

r« c l{Any^/'tld^x\g\'/'{Ti-D/2){rnlf/' 

+ (i - e)r(l - D/2){ml)^/^-'R} . (2.41) 



If we write 



r(i) = 5^rw 

n 

= y"d^xM^/2{4^+4^ + ---} (2.42) 

we have 

4^ = ^(4vr)-^/2rr(-Z?/2)^(m^)^/2, (2.43) 

n 

4) = \{4n)-^/'t{l-^)RTil-D/2)J2{ml)^/'-\ (2.44) 



The term in represents a contribution to the cosmological constant in 
D dimensions induced by quantum corrections to the classical theory. C^j^^ 
is recognized as the negative of the Casimir, or vacuum energy density, and 
corresponds to what we calculated earlier |23| in the case ^ = and D = 4. 



The term in >C^^ represents an induced gravity term which can renormalize 
the Newtonian gravitational constant found in the classical gravitational 
action. In the special case ^ = 1/6, which corresponds to the conformal 
value for the coupling constant when D = 3, we have = 0. In this case 
the lowest order correction to the classical gravitational action would result 
in an theory. For the original 5-dimensional Randall-Sundrum metric, 
D = 4, and the conformal value is ^ = 3/16, so there is an induced gravity 
term in the effective action. 



3 Renormalization 

The expressions for and ^JJ^ obtained in (|1|) and ( p4| ) may be 
evaluated using the method described in Appendix A or B. With uj{s) defined 
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in ( [A.ip we have 

/:« = h4Tr)-^/H'{ka)^uj{-D/2) 



(3.1) 



^(1) 



i(4^)-^/2r Q - i?(A;a)^-2^(l - D/2) . (3.2) 



The pole parts of these expressions wih determine the form of the required 
counterterms. If we use PP{- • •} to denote the pole part of any expression, 
from ( |A.45D and (|A.46| ) we find 

PP{4)} = 1(4.)-^/^^ 



T{D/2) 



a^ + (-l)^ 



d 



D 



(3.3) 



{D - 2)A;^-2 

r(i?/2) 



(3.4) 



The first few coefficients are given in (|A.27HA.30| ) for untwisted fields, 



and in (B.5-B.10) for twisted fields. In both cases higher order coefficients 
may be calculated in a routine manner, but because they become lengthy 
we will not give the results. 

The renormalization can be effected by modifying the procedure used by 



Ref. 1 24], which used the brane tensions. Because we allow curvature on the 
branes, it is necessary to augment the brane tension terms with curvature 
or else it will not be possible to remove the pole terms coming from (|3.4| ) 
when D is even. As noted in R ef. [p4| , the key feature is the a-dependence 
of the results. When D is odd, (|3.3| ) and (3^) show that the pole terms are 
proportional to — 1 and a^~'^ — 1 respectively. Such pole terms can be 
dealt with by a renormalization of the bulk gravitational and cosmological 
terms. The sign difference between the factors of a and the 1 comes about 
from the integration over the extra y coordinate to obtain the effective 4- 
dimensional theory. When D is even (D = 4 for the Randall-Sundrum 
model) the pole terms are proportional to + 1 and a^^"^ + 1 respectively. 
It is not possible to use the bulk terms to remove these divergences. The 
situation is easily understood from the viewpoint of heat-kernel methods 
and dimensional regularization, and was pointed out a long time ago ||28t] . 
In an odd spacetime dimension (as occurs for D = 4) the only possible 
counterterms come from the boundary of the spacetime. For the Randall- 
Sundrum model this corresponds to the two branes. 
The metric on the visible brane is 



9%{x) = g^uix,y = nr) 



(3.5) 



10 



The metric on the hidden brane is 

a^^Lvix) = 9iiu{x,y = (^) = g^,v ■ (3.6) 

These give rise to 

i?„ = a-'^R, (3.7) 
Rh = R (3.8) 

as the curvature scalars on the visible and hidden branes We will take 

Sbrane — ^brane ^brane i"^-^) 

where 

Slane = j d'^ x\g,\^'\K,R, - V,) , (3.10) 

Stane = j d"" x\gh\^'\KhRh - Vh) , (3.11) 

as the brane action, k^, Kh,Vy,Vh are all bare quantities. It follows that 

Sbrane = J d"" x\g\^/^ {(kh + a'"'^ K,)R - {Vh + o^K)} • (3.12) 

All bare quantities can be expressed in terms of renormalized ones and 
counterterms in the usual way : 

l'^v,h = + ^'^v,h 1 (3.13) 

Vv^h = + (3.14) 
This gives a counterterm brane action 6Sbrane from (3.12) in an obvious way. 



Because the structure of the pole terms in ( p. 3D and (3^) is the same as that 
of 6 Sbrane, there is no problem absorbing them into the counterterms. 

There is always the freedom to perform finite renormalizations since the 
renormalization counterterms should be fixed by imposing some renormal- 



ization conditions on the effective action. (See Ref. [29| for example.) To 
this end we can write 

SKv,h = —^ + ^i^v,h^ (3-15) 

SV,,h = -^ + 5V,\, (3.16) 
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where h^^^v h ^^^vh independent of e. Noting that the 

factors of in (3.12) are really a^'^'', we find 



{6kI + a^-25K° + a^~^6K-' In a) 



R 



rl 



+a^(5K° + a^'5Kr^ln( 



(3.17) 



as the Lagrangian density counterterms on the branes. D is now fixed to 
be the physical spacetime dimension (eg 4). In order that the poles in 
and (13.41) be removed we must fix 



6V' 



T{D/2) ' 

2)A:^-2 



T{D/2) 



5k 



-1 



(-1) 



1)^5kI 



(3.18) 
(3.19) 
(3.20) 
(3.21) 



If we use the results contained in the appendices for either twisted or un- 
twisted fields, (see ( |0^ , [05D ), it can be seen that involves either 
terms independent of a, proportional to or In a, apart from the term 
we have defined as Ai. Apart from Ai, the a-dependence is the same form 
as the finite counterterms SVj^ + a^6V^ + a^5V~' In a appearing in 6£brane- 
This means that all terms in apart from Ai can be absorbed by finite 
counterterms. A similar conclusion holds for the curvature dependent part 
of C. The renormalized expression for C\ becomes 



^Dy^R + i(47r)-^/2(A:a)^ A^{-D/2) 



(3.22) 



The renormalized expression for Cr becomes 

-(- 

Stt V6 



Cr 



K^+a^-2^^+-^(^-e)(4^)i-^/2(fca)^-2 Ai(l-Z?/2) R. (3.23) 



The function Ai{s) was defined in ( |A.37| ) with ( |A.21D used for gu{z) in the 
untwisted field case, and (B^) used for twisted fields. 

In the next section we look at the possibility of finding self-consistent 
solutions arising from quantum corrections to the classical theory. 
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4 Self-consistent solutions? 



Having obtained the vacuum, or Casimir, energy for the effective D-dimen- 
sional theory we can now see if there are any solutions of the Randall- 
Sundrum type to the quantum corrected field equations. We will argue that 
the problem is not quite as simple as calling the Casimir energy a potential 
and looking for its extremum. Rather, we must require that the Einstein 
field equations hold, and we will see that this has a different requirement. 

If r denotes the full effective action with quantum effects included we 
must require 

= 0, (4.1) 



5T 



: 0, (4.2) 

da 

with the Randall-Sundrum spacetime ( p.ip as a solution with g^^ = r]^u. The 
classical part of T comes from ( ^.4]) and (|2.5| ) , and to these contributions we 
must add the one- loop correction T^^\ If we perform the integration over 



the extra dimension y using the metric ansatz (2A) with a given by (|2^ 
we obtain 

Sg = (WttG)-' j dv,[ ^ (1 - a^-^)fi + AkD{l - a^) 

-2fc(I) + l)(l-a^)-i^(l-a^)}. (4.3) 

dvx = d^x\g\^^'^ is the D-dimensional invariant volume element. Because 
we set g^^i, = r}^^ after differentiation with respect to g^j_u, the term in the 
curvature R will make no contribution to the field equations. The same 
will be true for the induced gravity part of the one-loop effective action, 
as well as for higher order curvature terms that we have not calculated. 
The induced gravity term is important for relating the overall coefficient of 
R in the effective action to the physical value of the gravitational constant, 
which is irrelevant to the field equations in the present case. We are therefore 
justified in only considering the vacuum energy in the present calculation. 
Accordingly we will take 

T = Sg- j dvx F{a) (4.4) 

where F{a) includes the brane tension terms from ( p.Sp as well as the Casimir 
energy. 
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Variation of T with respect to g^jy followed by taking y^i, = r/^,y results 

in 

OA 

k{D - 1)(1 - a^) - —(1 - a^) = 87rGF{a) . (4.5) 
kD 

Variation of T with respect to a yields 

k{D-l)-^ = -'-^a^-^F'ia), (4.6) 



after setting g^^^, = rj^y. Combining ( |4.5| ) and ( [4.6D gives us a requirement 
for F(a) : 

Q = {l-a'^)aF'{a)+Da^F{a) . (4.7) 
As a check on this last result, if we set 

F{a) = Ka^ + Vh , (4.8) 

corresponding to the brane tension terms in ( |2.5| ), we find the condition 

V, + Vh = Q. (4.9) 

This reproduces the result on the balancing of the two brane tensions needed 
in the Randall-Sundrum model. The remaining requirements ( |2.7] ) and ( |2.8| ) 
follow upon using ( ^l9| ) in ( |45| ) and (^^. 

If we now include the vacuum energy by taking 

F(a) = Ka^ + + /(a) , (4.10) 

in place of the classical expression ( |4.8| ) , we must have 

= Da^{V^ + Vh) + (1 - a^)af'{a) + Da^f{a) (4.11) 

holding. This is clearly not just a simple minimization condition on the 
potential. From ( 3.22| ) we have 

/(a) = -i(4vr)-^/2A;^a^yli(-Z)/2) (4.12) 

with Ai defined in ( [A. 37 ). If we take 

I{a)=Ai{-D/2), (4.13) 

then we must have 

2L>(47r)^/2 , , , n^ ^ . . 

= ^-j^ (K + Vh) - DI{a) - a(l - a^)/'(a) . (4.14) 
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This condition must be fulfilled if the Randall-Sundrum spacetime is to be 
a self-consistent solution to the effective 4-dimensional Einstein equations 
with quantum effects included. 

In order to satisfy the self-consistency relation ( 4.14 ), we can use the 
freedom to tune the brane tensions. A possible choice is to require that the 
balancing condition ( [4. 91) between the brane tensions still holds. For this 
choice it is easy to see that there are no solutions to ( 4.14 ). To see this is 
quite simple. It is easy to show that /(a) is negative for < a < 1, and is 
monotonically decreasing. The right hand side of ( 4.14 ) can never vanish. 
However, there is no need to require that V^ + Vh = when quantum effects 
are included^, and relaxing this condition allows solutions to ( 4.14 ). The 
interesting thing is that it is possible to find self-consistent solutions for 
values of a which solve the hierarchy problem only at the price of severe 
fine tuning of the brane tensions^. This is a more general result than [^] 
who only were interested in self-consistent solutions which had values for a 
consistent with the solution to the gauge hierarchy problem (ie a « 1), 
and who only examined the Casimir, or vacuum, energy. 



5 Massless, conformally coupled fields 

The case of massless, conformally coupled scalars deserves special mention. 
Taking m? = and ^ = {D — 1)/{AD) in the order of the Bessel functions 
determined in ( 2.181) gives u = 1/2. Because Bessel functions of half-integer 



order can be given in terms of elementary functions, the analysis becomes 
much simpler than for general v. If we examine the twisted fields first, and 
use 

Ji/2{z) = \ — sinz , (5.1) 
Yii2iz) = -J—cosz, (5.2) 

' V TTZ 

it is easy to see that Fi/2{z) defined in ( |2.32| ) is given by 

~ , , 4sin(l — a)z 

Fy^iz) = ; ' . (5.3) 

^We would like to thank J. Garriga for having pointed this out to us. 
^This has been noted in for the case of massless minimally coupled and conformally 
coupled fields 
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We therefore have the mass eigenvalues given by 



irka 
(1^ 



n , for n = 1, 2, 3, . . . . 



We can evaluate the sum over modes in £j^'' directly using the Riemann 
C-function, without using the contour integral method of Appendix B. It is 
easy to see that 



(5.4) 



(1) 



:(4vr) 



-D/2 



Trka 



D 



T{-D/2)a-D) 



2^ 



(4vr) 



-D/2 



ka 



1 



D 



(5.5) 

T[^^^\C{D + l) . (5.6) 



In the last line we have let e ^ 0, since the result is analytic at e = 0, and 
used some of the relations of the T- and (^-functions to simplify the result. 
Although no infinte renormalizations are required, it is still necessary to do 
finite ones. 

As a check on this result, we can now apply the method of Appendix B. 
We need to know Ai{—D/2) which is given by 



Ai{-D/2) 



T{l + D/2] 



oo 



T{D/2) 



dz z^-^ In 



1 _ e-2{i-a)^ 



1 



-2z 



(5.7) 
(5.8) 



(In the second line we have integrated by parts and made use of the Bessel 
functions (5.1.51^) for u = \/2. ) The integral is easily evaluated in terms 
of the Riemann (^-function with the result 



Aii-DII) 



TT I 2 



C(D + !)[(!- a 



11 . 



(5.9) 



Note that when expanded in powers of a, Ai begins with order a as found 
in the case of general u. The result in (|5.9| ) gives rise to the same result 
as found by a direct summation of the mass eigenvalues, after imposing the 
same renormalization conditions, and provides a useful check on the method. 
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The untwisted conformally coupled scalar field case was also dealt with 
in Ref. |]2^ who exploited the conformal invariance to relate the problem 
to one in flat spacetime with boundaries present. It is worth emphasizing 
that this requires the theory to have no conformal anomaly, or else extra 
terms will arise when transforming back to the original spacetime. The 
fact that no such terms do arise is guaranteed by our method which makes 
no use of conformal transformations. In effect we have shown that there 
is no conformal anomaly for conformally coupled scalars in the Randall- 
Sundrum spacetime in any dimension. Vanishing of the conformal anomaly 
is guaranteed by the vanishing of when v = 1/2. (This is easily checked 
to be true for the coefficients listed in (B.5- BT0| ).) 



We now look at the case of untwisted fields. The functions j,y and Ui, 
defined in ( 2.25 , 2^26|) are easily seen to be 



ji/2{z) = cosz , (5.10) 

yi/2{z) = \ — sinz , (5.11) 



We then find 



^1/2(2) = -^2sin(l-a)z, (5.12) 

where Fi/2iz) was defined in ( |2.24| ). The mass eigenvalues are given by the 
same expression ( ^.41) as we had in the twisted field case. A direct summation 
over the mass eigenvalues again leads to ( |5.6| ). 

We can also check the result by using the contour integral method of 
Appendix A. The integral Ai{—D/2) becomes 

1 ^fD + l\ _o , ^ 



A,{-D/2) = — r ( ^- ) C{D + 1)[(1 - + 1 - 2-^] . (5.13) 



The only difference with the twisted field result is the constant term. The 
(1 — a)~^ dependence is the same. Note that as a consequence, when a is 
small, Ai is approximately a constant, rather than proportional to a. This 
would seem to contradict the general results established in Appendix A. 
The root cause of this lies in the behaviour of Fi/2{z) as 2; ^ 0. Clearly 
^1/2 (-2) ~ z^ as z ^ 0, which does not satisfy the assumption we made in 
Appendix A that -^1/2(0) is finite and non-zero. In this particular case we 
must remove the zero of Fi/2{z) at 2; = and take instead z~'^Fi/2{z) which 
does satisfy the conditions of being finite and non-zero as z ^ allowing 
the deformation of contour used in Appendix A. If we compare z^'^Fi/2{z) 
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with Fi/2iz) it is obvious that the same result wih now be obtained for 
Ai{—D/2) for the untwisted field as we found for the twisted field, and 
which is in agreement with the direct summation over mass modes. In any 
case, the difference between the twisted and untwisted fields for Ai results 
in a term proportional to in the effective action which can be dealt with 
by renormalization. 

For either twisted or untwisted field configurations it is easy to show 
that the only solutions to (4.14) is a = for Vy + Vh = 0. When this last 
condition is relaxed, it is possible to find self-consistent solutions, but if we 
demand the hierarchy problem to be solved a severe fine tuning of the brane 
tensions is required. 



6 Discussion and conclusions 

In the present paper we have examined the one-loop radiative corrections to 
the effective action arising from the quantization of a massless scalar field 
non-minimally coupled to the curvature of the background spacetime pro- 
posed by Randall and Sundrum [^. Although it is not possible to obtain an 
exact result for general, curved versions of the Randall- Sundrum spacetime, 
it is possible to obtain an expansion in powers of the curvature. The vacuum 
(Casimir) energy is the first term in this expansion. The vacuum energy has 
received attention due to its possible role in stabilising the ra- 

dius of the extra dimension in a self-consistent way, as an alternative to the 



mechanism proposed in |20|. Because of the different conclusions and in- 
consistent results obtained in previous analyses, we have presented a careful 
and complete calculation here. Differences and comparisons of our results 



are made with Refs. |^, 23, |2J] at appropriate stages. 

One important clarification concerns the boundary conditions adopted 
for the scalar field. We emphasized that there were two inequivalent field 



configurations (conventionally referred to as twisted and twisted [26|). Ei- 
ther can be used, and in general the results will be different. 

Another clarification concerns the renormalisation of the theory. We 
presented a full calculation for both field configurations for the generally 
coupled massive scalar field with all relevant counterterms indicated. To 
do this it is necessary to perform a full analytic continuation of relevant 
integrals as outlined in Appendices A and B. As we emphasized, the same 
method can be used for both twisted and untwisted fields, but the precise 
details of the analytic continuation are slightly different due to the different 
asymptotic behaviour of the Bessel functions between the two field config- 
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urations. Ref. [^] did not perform the detailed analytic continuation, and 
really only analyzed the untwisted field configuration. Nevertheless, as we 
showed above, the essential conclusions of their paper concerning the renor- 
malisation are correct in either case. One reason for presenting the detailed 
calculation of the counterterms is that the pole term in dimensional regu- 
larization is related to the conformal anomaly for the scalar field. 

As noted in the older Kaluza-Klein models based on spaces with homo- 
geneous extra dimensions [21 1, a full analysis of self-consistency can only be 
achieved if the quantum correction to the gravitational part of the effective 
action is known. We presented the necessary analysis for this above. In 
order to renormalize the theory when the branes were curved, it was neces- 
sary to add on additional contributions to the brane tensions involving the 
curvature. This suggest that the brane sector of the model has to be viewed 
as a series, with the brane tensions as the first term in the expansion. 

Another clarification concerns the conditions for self-consistency of the 
Randall-Sundrum spacetime. With flat branes we showed that the induced 
gravity contribution was not important, in contrast to the situation in |21]. 
Additionally, it was argued that the requirement of self-consistency is not ob- 
tained simply by looking at the minimum of the Casimir energy. By looking 
at the quantum corrected Einstein equations for the effective 4-dimensional 
(or more generally D-dimensional) theory we obtained a consistency condi- 
tion which the radius of the extra dimension had to satisfy. We found no 
solutions to this if we demand Vy + Vh = 0, meaning that it is not possible 
to stabilize the radius by quantum effects coming from bulk scalar fields in 
this case. This conclusion can be altered by relaxing the balancing condition 
Vv + = 0, even if we demand that the gauge hierarchy problem also be 
solved, but only at the price of fine tuning the brane tensions. 

We also examined massless conformally coupled scalars fields as a limit- 
ing case of our general results. This provides a useful check on the results, 
since it is possible to perform the calculations in a different, and much sim- 
pler, way. It also provides a comparison with Ref. |22] where a conformal 
transformation was used to relate the calculation to one in a flat spacetime 
with boundaries. A curious feature of the massless conformally coupled case 
is that the imposition of twisted or untwisted boundary conditions makes no 
difference to the final result. As we mentioned above, our calculation pro- 
vides a proof that there is no conformal anomaly for the Randall-Sundrum 
spacetime in any number of dimensions. 

An interesting point concerns the relation between viewing the quantized 
scalar field as an effective 4-dimensional theory, and a 5-dimensional theory. 
From the 5-dimensional viewpoint, self-consistency must be satisfied from 
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the 5-dimensional Einstein field equations. To do this, it is necessary to 
evaluate not just the total energy, but rather the energy density. (More 
specifically we must know the vacuum expectation value of the stress-energy 
tensor.) Due to the boundary conditions and the nature of the spacetime, the 
vacuum energy density for the 5-dimensional scalar field will not be expected 
to be constant in general. It would be interesting to study this problem, 
although the presence of Bessel functions in the mode eigenfunctions would 
make this rather difficult. 

Final comments concern the whole question of using quantum effects to 
obtain a self-consistent result for the radius of the extra dimensions. The 
conclusions of Ref. |2j] were rather negative concerning the possibility of 
simultaneously solving the gauge hierarchy problem and self-consistency. It 
is still interesting to ask if it is possible to obtain self-consistent solutions 
at all, regardless of the size of the extra dimensions. All that we currently 
know about this is that it is possible to do so with bulk scalar fields, but, fine 
tuning of the brane tensions is required when we demand that the hierarchy 
problem is simultaneously solved. 

Whether or not a theory with more realistic matter content can do so 
is an open question. Finally, even if it is possible to find a self-consistent 
solution there still remains the problem of determining if it is stable. 
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A Evaluation of the sums for untwisted fields 

In this appendix we outline the evaluation of the sum over modes needed to 
calculate the vacuum energy and induced gravity terms. Begin by defining 

oo 

u{s)=r{s)Y,^n'' , (A.l) 
n=l 

where Xn is the ra*^ positive solution to Fy{xn) = with Fi,{z) defined in 
terms of Bessel functions in ( ^.24| ). For large n we know x„ ~ n so that 



(A.l) is defined initially for 3f?(s) > 1/2 where the sum converges. We 
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want to analytically continue u;(s) to s = —D 12 and 1 — Z)/2. The basic 
method for doing this is to convert the sum into a contour integral. A simple 
application of the residue theorem shows that 

^{s) = ^ [dz z-'^ ^ In F,{z) , (A.2) 
ZiTi Jq dz 

where C is any contour which encloses the positive zeros of Fjj{z) defined in 
( |2.24| ). A convenient choice is shown in Fig. 




Figure 1: This shows the contour C used in the evaluation of (A.2). The contour 
is traversed in the counterclockwise direction and the hmits r ^ and R —> oo are 
understood. There is a branch cut in the integrand that we choose to he along the 
negative real axis. 



This procedure is basically the same as that used in Refs. ^ for 
studying the Casimir effect in a spherical shell. Ref. [^0| uses an integration 
by parts to remove the derivative in (A.2), and it is this basic starting point 
that Refs. |22, 24 1 follow. It is more convenient to stick with ( [A .21 ) for our 
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purposes, although either method could be used. A difference between our 



approach and that of Ref. |31| is that we do not need to use the uniform 
asymptotic expansion of the Bessel functions, making our calculation slightly 
simpler. 

By using the behaviour of the Bessel functions for small z it can be shown 
that the contribution from the semicircle of radius r behaves like r^~^* for 
small r. This means that if we restrict ourselves to the strip ^ < < 1 
the contribution to the contour coming from this portion of the contour 
will vanish as r — > 0. (The fact that -Fi/(0) is finite and non-zero is crucial 
to this result.) For the semicircle of radius ii, the integrand behaves like 
R~^^ for large i?, which follows from using the large argument expansion 
for the Bessel functions in Fyi^z). The contribution from this portion of the 
contour also vanishes if we let -R — > oo and restrict ourselves to the strip 
\ < 5R(s) < 1. The only contributions to the integral for ^ < 3f?(s) < 1 come 
from the portions of the contour along the positive and negative imaginary 
axis. After a short calculation it can be shown that 



U!{s) 



oo 



with 



and 



2 

P„(z) = - [iy{az)ki,{z) - iy{z)ki,{az)] . (A. 4) 



iu{z) = j{l-40l.{z) + zlUz) , (A.5) 
kuiz) = ^{l-Ai)Ky{z) + zKl{z) . (A.6) 

We now wish to analytically continue the expression (|A.3| ) for uj{s) out 
of the strip ^ < 5R(s) < 1 to 5R(s) < ^. The simplest way to do this is to 
note that the impediment to convergence of ( [A.3| ) for 3f?(s) < ^ comes from 
the behaviour of Pv{z) at large z. We will define 

Iy{z) = (A.7) 



Ku{z) = J-e-^s(^)(z), (A. 



which gives 




i.iz) = J-e^n"{z), (A.9) 
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ku{z) = -y^e-^sW(z), (A.IO) 
from ([A .51) and ([A.6|) where 



e(^)(z) = s(^)(z) + ^S(^)(z) + ^[I)(l-40-l]s(^)(-), (A.ll) 
= s(^)(z)-^sW(z)-^[Z)(l-40-l]S(^)(^)(.A.12) 
For large z the asymptotic expansions for the Bessel functions show that 

oo 

S«(z):^^a,z-\ (A.13) 

A;=0 

where 

_ (-i)fer(z. + fc+i) 

and 

~ . (A.15) 

Because the divergent behaviour of oo{s) comes from the large z be- 
haviour of the integrand, we will split up 

uj{s) = UJi{s) + UJ2{s) , (A.16) 

where 

1 /■ - -2s d 







oo 

^2{s) = dz z~^^ P,{z) . (A.18) 



(The choice of 1 as the upper limit in ( A.17 ) is not important; any value 
greater than will be equally good.) Analysis of the terms in Pu{z) respon- 
sible for divergence when 3f?(s) < ^ leads us to define 

P,{z) = -a^l^ze'^^--^ T.^J\z)T.\f\az)g,{z) , (A.19) 

TT 
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where 



9u{z) 



1 - e 
1 



-2{l-a)z '>^i^\z)^u\az) 

li!\z)T,f\az) 
ki,{z)iy{az) 
iu{z)ku{az) 



(A.20) 
(A.21) 



Substitution of ( A.19| ) into (A.l?) results in 



1 I- a 

+ 



r(l - s) i 2s 1 - 2s J r(i - s) 



+ 



1 



dz z 



+ ^ln 

dz 



]«(z)sf)(az)]} 



(A.22) 



where the result has been analytically continued to 5R(s) < 0. Substitution 
of ( [A.19| ) into ( [A.18D , assuming initially ^ < 5R(s) < 1, and then analytically 
continuing to K(s) < 0, results in 



w(s) 



r 1-s 



dz z '^'^ — In Quiz) 
dz 



+ 



1 



r(i-s) 
1 



dzz ^''-^In 
dz 



dzz In 



d_ 

dz 



S(^)(z)sf)(az) 



)(az)l . (A.23) 



The first two terms are analytic for 5R(s) < 0; however, the last term has 
poles at various values of s which we need to determine if we are to perform 
the required analytic continuation. 

We have constructed T,u (z) so that T.P{z) ~ 1 + O {z ^) as z ^ oo. 
This means that we can write 



In J:^J\z)^J2dk 



(A.24) 



k=l 



for large z, with some coefficients dk- The coefficients dk will determine the 
poles of Lj{s). Because of ( [A.15 ) we have 



In 



oo 



(A.25) 



k=l 
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with 



1 + 



(-1)* 



(A.26) 



Using ( |A.13 ) and (A. 14) in ( [A. 24 ) it is easy to calculate the coefficients dk- 
The first four coefficients are found to be 



di 
d2 
d3 
di 



3 1 2 

a u 

8 2 



—V 1 — a a 

4 16 2 2 

1 2 3 



23 



-—av 
2 



-a" + -a H V H v 

2 8 48 24 



21 1 . 

h -a^ 

128 3 



1 

— ( 
2 



1 2 3 



13 



—a H — a H 

2 8 16 



where 



-a 



4 
D 



27 
128 



(l-4e). 



(A.27) 
(A.28) 
(A.29) 

(A.30) 
(A.31) 



(For a = 2, which is the case for D = 4 and = 0, the results reduce to 
those given in (2^.) By adding and subtracting terms to the last integrand 
of ( A.23| ) we may analytically continue the result to < 0. The more 
negative 3f?(s) becomes, the more terms we need to subtract. In this way we 
obtain 



n 



oc 

1^/ 



dz z 2« -^^ In gy{z) 
dz 



+ 



+ 



ri 



dzz ^''-^In 
dz 



1 



r(i-.) 
1 



dzz-"^'— I In 



dz 



^i'\z)^f\az 
S(^)(z)sf)(az 



N 



kcTk 



N 

E 

k=l 



akz 



-k 



(A.32) 



Provided that N is chosen so that ?R.{N+l+2s) > the analytic continuation 
to any value of s with 5R(s) < may be performed. All of the possible poles 
can be found in the last term of ( A. 321) . Some procedure like the one we 
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have described is essential to provide the proper analytic continuation of 
cj(s) to < 0. 

We are interested in the value of a;(s) at s = p — D/2 for p = 0,1. There 
is a simple pole at this value of s coming from the k = D — 2p term in the 
sum in the last term of ( [A.32| ). We can write 

uj{s) = u;poio(s) + uJrcgis) (A.33) 

where 

"-^^^'^ = " r(i-.) (z;-2^ + 2.) 

contains the pole term (as well as a finite term found by expanding about 
the pole), and a;reg(s) is regular [ie analytic) at s = p — D/2 and is given 
by an expression like ( [A. 32 ) but with the term k = D — 2p omitted from 
the sum in the last term. We are free to simply set s = p — D/2 in ujj.Qg{s), 
where D is the physical dimension (ie 4 for the Randall- Sundrum model). 

The most important thing about the regular part ujj.(,g{s) is the depen- 
dence on a. By splitting up the logarithms in ( |A.32D it can be seen that 



Wreg(s) = C{s) + Ai{s) + A2{s) (A.35) 



where 



1 



-j^ \ k — 1 



r(i-5)^ ik + 2s) 

is independent of a, and 



oo 



1 r ^ ~2s d 





1 
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7 ( N 



4(-i)* 



contain the a-dependence. The prime on the sums in ( A.36| ) and ( A.38| ) 
denote that the term k = D — 2p, which gives rise to the pole, is omitted. 

The a-dependence of ( [A.38| ) can be made more exphcit by a change of 
variable z ^ z in both integrals : 





oo 

,2s 



k=l 



N 



oo 1 oo 

Now split up the integration range of the second integral as / = / + 



a a 1 

After a short calculation we are left with 



r(l-s)^ (fc + 2s)a*^ ^ ' 



where 
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is independent of a (and is analytic for < 0). Finally we may perform 
the integration in ([A .401) to find that 

A2{s) = a'^A{s) - ^-^M_^^dD-2pa''-''lna (A.42) 



with 



We can conclude that for s = p — D /2 

u;regip-D/2) = C{p-D/2)+Ai{p-D/2)+a^P-^A{p-D/2) 

T(l-s)^^ [k + 2s) ^ ' 



where C(s), j4(s) are independent of a. 

Now turn to a similar analysis for the simpler pole part of given in 
( [A.34 ). Because there is a pole, we will use dimensional regularisation and 



set D — > Z)+e, with s = p—{D+e)/2. Note that it is only in the s-dependent 
terms that we replace D D+e, since the factor of {D—2p)a£)~2p came from 
a term in the asymptotic expansion which involved only integral powers. 
Because ak oc dk, and depends on the order of the Bessel functions 
z^, there can be a D-dependence hidden inside aD-2p- However it is easy 
to see that this dependence can only affect terms which are constant, or 
proportional to a?P~^ . Such terms can be absorbed into redefined functions 
C{p — D/2) and A[p — D /2,) respectively, which we show in Sec. ^ can be 
removed by finite renormalizations. There is no point in giving the explicit 
form of such terms, and we can safely ignore the e-dependence of aD-2p- It 
is possible to write 

^poie(s) = ^ + ao + ai a^^'^ (A.45) 

for 



with ao and ai independent of a. (As explained above, it is a straightforward 
matter to obtain explicit results for qq and ai at s = p — D/2 hy expanding 
( [A^ ) about the pole.) 
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B Evaluation of the sums for twisted fields 



In this appendix we look at the calculation of w(s) using the solution to 
( |2.31| ). The details are very similar to those of the preceding appendix, so 



we will be brief. uj{s) is expressed as a contour integral as in ( A. 2 ) where 
F^{z) appears in place of F^{z). ( [A.3| ) still holds, but with Pu{z) now given 
by 

Puiz) = - [Iu{z)K^{az) - h{az)K^{z)] . (B.l) 
vr 

Using (^) and (^) results in 

P^{z) = -a-^/^z-^e^^^''^^T.'~^\z)T.'~^\az)g^{z) (B.2) 



vr 



where 



Ky{z)Iy{az) 

9v{z) = 1 - , . , . ■ (B.3) 
K^{az)I,y{z) 

The difference between the twisted and untwisted cases lies in the asymptotic 
behaviour of the integrand of the contour integral for large \z\. (Compare 
( [Al9| ) with (p^.) In the untwisted case the leading term in the asymptotic 
expansion comes from the derivative of the Bessel functions in ( |A.5 , A.(^ ). 



These terms are absent from the twisted field expressions. Because of this 
different asymptotic behaviour, which lead to divergent integrals dealt with 
by adding and subtracting terms as discussed in the previous appendix, 
it is difficult to deal with both twisted and untwisted fields in one single 
calculation. 

The steps in the analytic continuation of uj{s) follow much as before 
with slightly altered expressions. In place of ([A.23| ) we have the simpler 
expansion 

oo 

lnS(^)(z) ~ ^4^"'' (B.4) 
k=l 

where the first few coefficients are given by 

di = ^(1-4^.2), (B.5) 

d2 = 3^(1- 4^') > (B.6) 
ds = ^(25 - 104Z.2 + 161.4) , (B.7) 
dA = Y^(13-56i/2 + l6j.4) , (B.8) 
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4 = ^^(1073 - 47481/2 + I840i/^ - QAu^) , (B.9) 
5120^ ^ ^ ^ 

(Je = Y^(103 -465i/2 + 216i/^-16z/^) . (B.IO) 

The remaining analysis of Appendix A follows as before with these new 
coefficients in place of the dk, and ( p. 3D used for Quiz). (This changes the 
result we called Ai(s) in ([AIstI ).) 
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